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Abstract 



We continue our examination of the constraints in spherically symmet- 
ric general relativity. We extend to general configurations with J ^ 
the analysis of II which treated a moment of time symmetry. We exploit 
the one parameter family of foliations introduced in I which are linear 
and homogeneous in the extrinsic curvature to characterize apparent 
horizons and spatial singularities in the initial data. In particular, we 
demonstrate that these characterizations do not depend sensitively on 
the foliation. 
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1. INTRODUCTION 



This is the third paper in a series in which we examine the general features of the con- 
straints in general relativity under the assumption that the spatial geometry is spherically 
symmetric and possesses just one asymptotically flat region [1,2]. In paper II, we focused 
on solutions of the constraints which occur when the extrinsic curvature Kah momentarily 
vanishes (MSCs). As such, we did not need to address the issue of fixing the foliation. 
In this paper, we extend this analysis to incorporate a non-vanishing extrinsic curvature. 
Unless one sets out to be difficult, this corresponds to a non-vanishing flow of matter, J. 

The introduction of extrinsic curvature complicates the analysis substantially. The 
advantage of having dealt separately with moment of time symmetry configurations in 
paper II is that we can focus here on the physical feedback on the spatial geometry intro- 
duced by extrinsic curvature. The important point is that the solutions of the constraints, 
as well the relationships between the global measures of the energy and the dimensions of 
the support of matter which we exploit to characterize horizons and singularities, are not 
sensitively dependent on the gauge fixing the foliation. 

To fix this foliation we implement explicitly one of the gauges parametrized by a, 
linear and homogeneous in the extrinsic curvature, introduced in paper I. These gauges 
will serve to set in context our understanding of the constraints when the initial data is 
momentarily static [2]. It was shown in paper I that the allowed values assumed by the 
parameter correspond to all tangent vectors lying within the superspace lightcone. One 
of these gauges is the maximal slicing gauge. Another is the polar gauge. Among the 
attractive properties of all such gauges is that they foliate flat spacetime by flat spatial 
hypersurfaces. When the momentum constraint is satisfied, the extrinsic curvature is linear 
in J, albeit in a non-local way. In this way the extrinsic curvature of the hypersurface 
responds directly to the (radial) movement of matter on it. Another possibility, which is 
motivated by the introduction of the optical scalars as canonical variables on the phase 
space, is to treat the trace of the extrinsic curvature, itself as an independent datum 
along with energy density, p, and J. While this is a legitimate gauge, it is not a usual one 
unless K — for ii K ^ the extrinsic curvature cannot adjust itself to the movement of 
matter. For this reason it does not correspond to our physical expectations and, therefore, 
we do not consider this possibility further here. 

As we did in paper II, we will focus again on the identification of the global structures 
that characterize non-trivial geometries — apparent horizons and singularities. 
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The identification of apparent horizons is comphcated by the fact that these physical 
landmarks no longer coincide with the extremal surfaces of the spatial geometry as they 
do at a moment of time symmetry: the movement of matter generates extrinsic curva- 
ture, thereby affecting how the spatial geometry is embedded in spacetime which, in turn, 
determines the lightcone structure on the surface. 

We begin in sect. 2 with a discussion of the generic analytic structure of the constraints. 
We derive a spatial diffeomorphism invariant analytic expression for the behavior of the 
geometry in the neighborhood of a generic singularity. Generally, the singularities of the 
three-geometry consistent with the constraints will be more severe than those which are 
possible at a moment of time symmetry. If, however, the movement of matter is tuned so 
that the extrinsic curvature vanishes as the singularity is approached, the strength of the 
singularity will be determined entirely by the quasi-local mass (QLM), exactly as it was 
at a moment of time symmetry [3]. This tuning corresponds to an integrability condition 
on the current. If, in addition, the tuning is refined such that the QLM also vanishes as 
we approach the singularity the curvature singularity disappears and the spatial geometry 
pinches off in a regular way. This latter integrability condition involving the QLM is 
completely analogous to the integrability condition we encountered at a moment of time 
symmetry. Regularity at the singularity is, of course, precisely the condition that the 
interior be a regular closed universe. If the matter fields carry conserved charges these 
will, in their turn, have integrability conditions associated with them. Viewed this way, 
regular closed universes appear to be very special universes [4]. 

In paper I, we represented the configuration space of the spherically symmetric theory 
by bounded closed trajectories on the optical scalar plane. In sect. 3 we examine these 
trajectories in vacuum. We discover that any trajectory that finds itself outside a proper 
subset of this domain is necessarily singular. While the detailed structure of this good 
subset depends on the specific gauge choice we make to determine the slicing, we find 
many features that are independent of the slicing. For example, on the boundary of the 
good subset we find two unstable fixed points. These correspond to the situation where the 
exterior spatial geometry neither collapses to a singularity nor expands to be asymptotically 
flat. Instead it becomes a semi-inflnite cylinder of fixed radius. 

We follow paper II by establishing global necessary and sufficient conditions for the 
occurence of apparent horizons and singularities. These conditions are framed in terms 
of inequalities which relate some appropriate measure of the material energy content on 
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a given support to a measure of its volume. The challenge is to identify useful measures 
in both cases. In analogy to the total material energy M (defined in paper I), we can 
introduce the total momentum, P, corresponding to the integrated material current over 
the proper spatial volume. The sufficiency criteria for the formation of a future (past) 
apparent horizon can be cast in a form which is a straightforward generalization of the 
moment of time symmetry inequality: if the difference (sum) of the material energy and 
the material momentum exceeds some universal constant times the proper radius, £o, of the 
distribution, the geometry will possess a future (past) apparent horizon. The corresponding 
constant for singularities is larger but the inequality docs not involve P. As wc found at a 
MSG a more appropriate measure of the material energy for casting the necessary criteria 
is the maximum value of the energy density of matter, pMax- The obvious generalization is 
the sum pMax + -^Max- However, we find that that the inequality is not symmetrical under 
interchange of p and J. If the dominant energy condition is satisfied, however, we can 
however cast the inequalities in the momentarily static form: if (pMax + <^Max)^o < some 
constant, the distribution of matter does will not possess a singularity for one constant 
and an apparent horizon for some other smaller constant. These inequalities are new. 

In the treatment by Bizon, Malec and 6 Murchadha (BMOM), and more recently 
by Malec and O Murchadha of the sufficiency conditions, the slicing of spacetime was 
always assumed to be maximal with K = [5,6]. If these inequalities are to be interpreted 
physically, they should, at least qualitatively, be reproduced in other gauges. We note 
that we never needed to address this question in paper II because the notion of a moment 
of time symmetry is gauge independent. We examine the sensitivity of these inequalities 
on the value of a appearing in the gauge condition. Not surprisingly, we find that the 
strength of the corresponding inequality does depend on the foliation gauge but not in any 
significant way so long as we are not close to the lightcone in superspace. 

Unfortunately, unlike in our examination of the constraints at a moment of time 
symmetry where we could fall back on the piecewise-constant density models, we enjoy no 
such exactly solvable standbys here. Even the analogue of the constant density star proves 
to be analytically intractable when J ^ 0. It is not surprising therefore that it is far more 
difficult to identify sharp inequalities than it was at a moment of time symmetry. Much 
of our effort is spent bootstrapping on moment of time symmetry inequalities. 

There is one extremely useful exactly solvable model consisting of a moving shell. We 
exploit this to speculate about the likely form of a possible generalization to J 7^ of the 
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lower bound on the binding energy derived by BMOM. 



2. THE CONSTRAINTS 

In this section we examine various general features of the constraints when Kab 7^ in 
a manner which parallels, wherever applicable, the treatment in paper II. We recall that 
the constraints are given by 

Kr [Kr + 2Kc] - [2 {RR')' - R'^ - 1] = Sttp (2.1a) 

and 

K'R+^{KR-Kc)=4nJ, (2.16) 
where the line element on the spatial geometry is parametrized 

ds^ = df + R^dQ^ , (2.1c) 

and we have expanded the extrinsic curvature (n" is the outward pointing unit normal to 
the two-sphere of fixed £) , 

Kab = nanbKc + {gab-nanb)KR. (2. Id) 

All derivatives are with respect to the proper radius of the spherical geometry, £. The 
spatial geometries we consider consist of a single asymptotically flat region with a regular 
center, £ = 0. The appropriate boundary condition on the metric at £ = is then 

R(0) = 0. (2.2) 

We recall that -R'(O) = 1 if the geometry is regular at this point. We assume that both p 
and J are appropriately bounded functions of £ on some compact support. A non-singular 
asymptotically fiat solution defined for all £ > will not, however, always exist for every 
specification of p and J. Our task is to understand what can go wrong. 

To solve the constraints classically, we need to implement a foliation gauge. This 
involves some spatial scalar function of the extrinsic curvature tensor. In a spherically 
symmetric geometry, this tensor has only two independent components. The foliation 
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either fixes one of these or relates it functionally to the second. Modulo the gauge, the 
momentum constraint can be solved for this other component. The extrinsic curvature is 
then completely determined by the sources. 

In paper I, we introduced the one-parameter family of gauges, defined for each a. by 



K£ + aKn = . (2.3) 

We showed that the momentum constraint can be solved uniquely in terms of the radial 
fiow of matter, J, as follows (Eq.(3.3) of I) 

KR^^J^d£R'+-J. (2.4) 

Whenever 0.5 < a < oo the gauge is valid everywhere and displays the correct asymptoti- 
cally flat falloff outside the support of J if the geometry is non-singular. When Eqs.(2.3) 
and (2.4) are substituted into Eq.(2.1), we obtain a second order singular non-linear integro- 
ODE for R* Subject to the boundary condition, (2.2), the solution is uniquely determined. 
Not only is the extrinsic curvature completely determined by the material sources, so also 
is the spatial geometry. There are no independent gravitational degrees of freedom. 

We note that in the gauge Eq.(2.3), the spatial geometry does not depend on the 
global sign of J. 

We saw in paper I that if Kr is regular at the origin then it must also vanish there. 
In fact, in the neighborhood of £ = 0, i? ~ £, so that 

Kur^AiT^^e. (2.5) 

To determine the n^^ derivative of at £ = we need to differentiate Eq.(2.1) n — 1 
times. A consequence of the vanishing of -fCit(O) is that J will only show up at order five 
— two orders behind p (see paper II). The behavior of the metric at the origin is clearly 
not sensitive to the current flowing there. 

It is instructive to also examine the values assumed by the optical scalars in the 
neighborhood of £ = 0. Recall that [1,7] 



* It is possible to rewrite Eqs.(2.1) so that they can be differentiated once to yield a 
Zoca/ third order singular ODE modulo Eqs.(2.3) and (2.4). However, it is not particularly 
illuminating to cast them this way. 
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UJ+ =2 R' + RKr , 



(2.6a, 6) 



LO- =2 R' - RKr , 



We can combine Eq.(2.4) of paper II and (2.5) to obtain 



a;±~2-^(p(0)T-^J(0)K^ 
6 2 + a 



(2.7) 



Suppose that the dominant energy condition (DEC) p > | J| is satisfied. If a > 1, then 
oj± < 2. If, however, a < 1 this is not the case. This demonstrates exphcitly that the 
inequahties (6.2a) and (6.2b) in paper I cannot generally be relaxed to the K = value. 
We note also that Eq.(2.7) implies 



which is consistent with the inequality (5.2) in paper I for all values of a. Note also that 
the absolute maximum of the product a;+a;_ obtains at the boundary values £ = and 
£ = oo and it is also the flat space value. When K = 0, this is also true of both oj^ and 
CO-. In general, the absolute maximum of neither need occur at these points. 

2.1 The Quasi-Local Mass 

As we found in paper II in a simpler context, the definition of the quasi- local mass 
can be exploited to provide an extremely useful first integral of the constraints. We recall 
that (Eq.(4.7) of paper I): 



Modulo the constraints (Eq.(4.8') of paper I), and the boundary condition, (2.2), 



If the geometry is non-singular, m is positive everywhere, coinciding in the limit ^ — > oo 
with the ADM mass, moo- Eqs.(2.8) and (2.9) are gauge invariant. To exploit Eqs.(2.8) 
and (2.9) to solve the constraints, we substitute the solution of the momentum constraint 
(2.4) in the gauge (2.3) into Eqs.(2.8) and (2.9). 





(2.8) 




(2.9) 
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Note that the leading spatial derivative in Eq.(2.8) is R' . Outside the sources, m is 
constant while ~ C/R}-^°'. From a functional point of view, Eq.(2.8) is identical to 
the energy integral in classical mechanics. To exploit this analogy, we recast Eq.(2.8) for 
all I as follows: 



where m is given by Eq.(2.9). 
2.2 The Neighborhood of Singularities 

In paper I, as a lemma to the positive quasi-local mass theorem, we proved that when 



if R > 1 anywhere the geometry must be singular. Let us suppose that R > 1 at some 
point. Then, when Kab satisfies Eq.(2.3) and a > 0.5, Eq.(2.1b) implies that R" < 0, so 
that R' is decreasing there. This can only occur by R' falling through R' = —1. Once 
R' falls below this value it will continue decreasing monotonically thereafter. The surface 
i?' = — 1 in the configuration space therefore acts as a oneway membrane. Suppose that 
the circumferential radius is Rq when R' = —1. We know now that the solution must 
crash, i.e. i? ^ in a finite proper distance which is less than or equal to Rq from that 
point. In fact, this is the only way the spatial geometry can become singular. 

How do we know that we have covered all possible singularities? We argue that the 
converse of the lemma holds. In general, — 1 < i?' < 1 if and only if the geometry is non- 
singular. As we will see, a singularity with i?' = — 1 is a result of a very special fine-tuning 
of the matter distribution. 

In the neighborhood of the point £ = £5 at which R = 0, Eq.(2.4) implies that 




(2.8') 



Now, formally at least. 




the weak energy condition p > holds, R^ < 1 everywhere in any regular geometry. Thus 



(2.10) 



Jll+a 



where 
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Ca(£) =47r / deJR^+'' (2.11) 
Jo 

is finite. Ku will therefore be singular (for any physically acceptable value of a) if the 
geometry pinches ofi^ unless the current is tuned such that 



C«(£s) = 0. (2.12) 

Now, if Kji 7^ 0, and a > 0.5, the most singular term in Eq.(2.8') is the quadratic in Kr. 
This implies that 

R'^ ~ R^Kl (2.13) 

in the neighborhood of R = 0, or R'^ ~ C^/i?^". Generically, therefore, R'^ diverges. The 
solution is 

f-%TV^' i^s-i)^. (2.14) 
\a + lj 

If q; > 0.5, such spatial singularities are more severe than the strong singularities discussed 
in paper II which are consistent with the Hamiltonian constraint at a moment of time 
symmetry. We will refer to the generic kind of singularity driven by extrinsic curvature 
as a strong J-type singularity. As a increases, the power law determining the strength of 
the singularity increases. Note that the limit a — > oo (the polar gauge discussed in I) is 
extremely singular. This is, however, a gauge artifact reflecting how poor the polar gauge 
really is. 

Unlike the strong singularities occurring in MSCs, at which the scalar curvature TZ 
remained finite, TZ will generally blow up (as Kj^ ~ l/(^s ~ ^)^)- dimensional grounds, 
we expect all curvature scalars to blow up as I /{is — i)"^ as we approach a singularity 
unless there is some constraint obstructing them from doing so. 

It is important to confirm that m remains suitably bounded as we approach a strong 
singularity. We do this by demonstrating that the volume integral (2.9) is always finite. 
We note that for suitably bounded p and J, 

(pR^r'.jr^Kr) ~ (p, j)(es-e)-i^) . (2.15) 
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If Q! < 2, the integrand itself remains finite. In general, the integral will be finite if the 
exponent of {Is — (■)~^ is bounded by one. But {a — 2)/(q! + 1) < 1 for all finite values of 
a thus guaranteeing that the integrals over E?pR' and R^JKn converge. 

It is clear that m{is) is always finite. Its sign, however, will depend on the details 
of the current flow. This is obvious from the definition Eq.(2.8). Even if R'^ > 1, a 
sufficiently large value of Kfi can render m positive. In particular, unlike the value of m 
assumed at strong p-singularities of MSCs which is always negative, the sign can assume 
either value. Indeed m need never even be negative in a singular geometry. Though 
R' decreases monotonically, m nonetheless remains positive. There is no conflict with 
the positive QLM theorem. In our examination of MSCs in paper II, we found that m 
is positive everywhere except at the origin or in a neighborhood of it if and only if the 
geometry is non-singular. This is a consequence of the coincidence of the converse of the 
bounded R'^ lemma and the converse of the positive QLM theorem when Kab = 0. In the 
general case, when Kab 0; such, coincidence occurs. 

What are the implications of the integrability condition, Eq.(2.12)? If Eq.(2.12) is 
satisfied the strong J singularity is moderated to one which is only strong a la p. The 
behavior in the vicinity of the singularity will then be determined by the m/R term in 
Eq.(2.8') even if the system was originally 'driven' towards the singularity by extrinsic 
curvature. If, in addition. 



the singularity will be a weak one with R'{£s) = —1- The corresponding bag of gold will 
be a regular closed universe. 

These integrability conditions depend on a. If a given function J satisfies Eq.(2.12) 
with one value of a, generally it will not satisfy that condition with any other value. What 
is missing is a spacetime diffeomorphism invariant statement of the integrability. 

If J is positive (or negative) everywhere, Ca(^) defined by Eq.(2.11) cannot vanish. 
Thus, if matter is collapsing or exploding everywhere, all singularities must be strong 
J- type singularities. 

This contrasts with the obstruction, p' < 0, discussed in paper II, prohibiting the 
formation of any singularity when Kab = 0. In general, we note that on performing an 
integration by parts on the first term, m can be rewritten 




(2.16) 
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m 








diR^ [JKr - p'] 



(2.17) 



The first term is manifestly positive. So is the third if p' < 0. If J is positive (negative) 
everywhere then so is m in any a - gauge. However, the third term appearing on the 
RHS of Eq.(2.8') may still pull the geometry into a singularity if J is sufficiently large. 
The peculiarity of momentarily static configurations with p' < discussed in paper II can 
clearly be destabilized by the motion of matter. 

All regular closed cosmologies simultaneously satisfy two integrability conditions, 
Eqs.(2.12) and (2.16). There can be no net flow of material from one pole to the other. In 
particular, J must change sign between the poles. In addition, Eq.(2.17) tells us that 



^0 

In particular, JKu — p' must change sign between the poles. These conditions will be 
examined in the closed cosmological context in a subsequent publication [4] . 

2.3 No strong J singularities in the Euclidean Theory 

The singularity structure we have investigated has one important consequence for Eu- 
clidean general relativity. If the sign of the quadratic term in Kji appearing in Eq.(2.8) 
had been negative, instead of facilitating the occurrance of singularities it would have 
presented an obstacle to their occurrance. Any non-vanishing extrinsic curvature would 
therefore tend to stabilize the spatial geometry against singularity formation. We note that 
there is precisely such a sign switch in the Hamiltonian constraint of Euclidean general rel- 
ativity. The Bianchi identities there tell us that the solutions of the constraints represent 
all possible configurations the system may assume as it is evolved with respect to Eu- 
clidean time. This suggests that gravitational instantons will tend to be more regular than 
their Lorenzian counterparts. In fact, the most singular Euclidean geometries will occur 
when the geometry is momentarily static. In a tunneling Euclidean four-geometry, such 
three-geometries correspond to the initial and final hypersurfaces of the Lorentzian space- 
times between which it interpolates. If these hypersurfaces are themselves non-singular, 
i.e. do not involve Planck scale structures, then Planck Scale physics does not enter the 
semi-classical description of tunneling between them. This would appear to validate the 
application of the semi-classical approximation. 




(2.18) 
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3 SOLUTIONS OF THE CONSTRAINTS AS TRAJECTORIES ON THE 
{uj+,oj+) PLANE 



In paper I we found that a very useful representation of the phase space was provided 
by the representation of solutions to the constraints as trajectories on the optical scalar 
plane. We have [1,7] 



(3.1a, 6) 



[iv+y = - 8ttR{p - J) - ^ f2w+ - 4 - 4UJ+KR - u+ujJ] , 
(a;_)' = - 87ri?(p + J) - 2^1-4: + 4lo_KR - a;+a;_] . 

3.1 Non-Singular Geometries off the support of matter 

In this section, we will focus on the behavior of trajectories outside the support of 
matter. We will suppose that the interior solution is regular. This is a more useful exercise 
than it might appear at first sight. This is because, as we have seen, the behavior of 
trajectories depends non-locally on the sources, p and J. In particular, the appearance of 
singularities does not depend sensitively on the values of p and J in the immediate vicinity 
of the singularity. In addition, as wc will see, the behavior of vacuum trajectories upon 
entering into a shell of matter is described in a simple way. 

Let us first recall briefiy the case the momentarily static solution outside matter. We 
set J = everywhere. We found that 

i?'^ = l-^. (3.1) 



where 



rrir 



Att / dm^pR' . (3.2) 



If moo < and R'(Iq) < 0, the vacuum geometry will be singular. If R'{£o) > 0, the 
vacuum geometry is non-singular but the positive quasi-local mass theorem tells us that 
the interior must harbor a geometrical singularity. 

If moo > 0, there is no way the geometry can be singular. If R'{£o) > 0, the solution 
grows monotonically. If R'{io) < 0, R will decrease until an apparent horizon forms at 
R = 2moo and thereafter increase. 
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When J the solution is considerably less simple. Let us rewrite Eq.(2.8) and (2.9) 

as 



R'^ = l- V{R) , (3.3) 
where the potential V{R) is given by 

V(R) = ^ - ^ , (3.4) 



rrir 



= 47r / MR^ [pR' + JRKr] , (3.5) 
Jo 



and Cot = Cai'^o) is given by Eq.(2.11). We suppose that 0.5 < a < oo. 

If moo is negative and R'{£o) < then, as before, the potential is monotonic and 
unbounded from below, and the geometry will be singular. It will generally be a strong 
J-singularity unless J is tuned such that Ca{io) = 0. If R'{io) > 0, the vacuum geometry 
is non-singular but the positive quasi-local mass theorem tells us that the interior must be 
singular. 

What is much more interesting is the case nioo > 0. Unlike the case of a MSG, a 
positive moo does not guarantee a non-singular exterior geometry. 

If moo is positive, the potential possesses a maximum at the point 

The value assumed by the potential at this point is 

V{R,) = {2-a-')^. (3.7) 

There are two possibilities we need to consider: 

If V{Rc) < 1, the geometry will be singular outside if and only if R'{io) < 0. 

If V{Rc) > 1, the nature of the geometry will depend not on R'{£o) but on the relative 
values of R{£o) and Rc- If R{£q) < Rc the geometry will always be singular. If R{£o) > Re, 
however, there is no way the geometry can be singular. The qualitative dependence on the 
sign of R'{£q) is the same as that for moo > when J = 0. 

The condition V{Rc) < 1 is equivalent to the inequality 
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< (2« - 1) ^^^J Cl. (3.8) 

Typically Eq.(3.8) will hold for a given moo whenever Cq, is large which corresponds, 
roughly speaking, to a large material current. However, because moo itself also involves 
J this criterion is not very precise. We need to isolate the dependence of moo on Kr. 
The inequality can then be cast as an inequality between R' and RKr. We recall that 
Ca = R'^RKn. Now Eq.(3.8) can be cast in the form 
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" (1 - R" + R'Kn^r < (2« - 1)^/' (^^) " l^^^l • 
We can also represent the inequality, R < Rc, as the exterior of the ellipse 



(3.9) 



R'^ + {2a - l)R^Kfi = 1 (3.10) 

on the (i?', RKn) plane. The beauty about Eqs.(3.9) and (3.10) is that when they are cast 
in terms of the optical scalar variables, a;+ and U- they are independent of R. 
Let us examine Eq.(3.9) in greater detail. We first cast it in the form 

f{RKR)<R'\ (3.11) 

where we define 



f{x) =x^- 2{2a - 1) (^^) + 1 . 

We note that f{x) is positive everywhere. In particular, /(O) = 1 and 



(3.12) 



/(±1/a/2Z^) = = f'{±l/s/2^r^) . (3.13) 
In the neighborhood of the points x = ±l/\/2a — 1, 



f^/'^ix) ~ \ '^^—^ \x T l/^/2^^\ . (3.14) 
\ a ' 

As a; — > oo, /^/^(a;) ~ x. The two branches of the function R' = ±f{RKjiy/^ correspond 
to the boundary V{Rc) = 1. They can be represented on the (a;+,a;_) plane (see fig. (3.1)) 
as the union of arc segments {CV, VP, PV', V'D} and {C'V, VQ, QV, VD'}, where the 
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coordinates of the points V and 7^', corresponding to the two mimima of /, are given repec- 
tively by (-2/V2a - 1, 2^2a - 1) and (2/V2a - 1, -2V2a- 1). Fig.(3.1) corresponds to 
a = 2. 

The eUipse defined by Eq.(3.10) is also represented on fig. (3.1) for a = 2. We note 
that for each a, the points V and V both he on this eUipse. The inequality, R < Rc, is 
represented by the region on the phase plane outside the ellipse. 

What is this figure telling us? There is a wedgelike region Qq, bounded by the arc 
segments, CQ, QD indicated on fig. (3.1) which determines the maximum excursion a 
vacuum trajectory can make from its point of departure, P = (2, 2), and still return home. 
This is a disjoint union of two regions, one in which V{Rc) > 1 and R > R^ the other in 
which V{Rc) < 1 and R' > 0. 

Any trajectory which lies outside Qq on exiting the support of matter is necessarily 
singular. This region, likewise, decomposed into a disjoint union, one in which V{Rc) < 1 
and R' < 0, the other in which V{Rc) > 1 and R < Rc- We note that these considerations 
did not rely on any energy condition, dominant, weak or otherwise. 

When the DEC holds, we note that Oq reduces to a proper subset of the domain, 
O, introduced in paper I to which all non-singular trajectories are confined For a = 2, O 
is given by the square, \u>±\ < 2. The region, — fio is rendered forbidden outside the 
support of matter. In particular, we note that the barriers ll!± = —2 are completely out 
of bounds. There always exists, however, a suitable p and J, which can be added within 
the region Q so as to render the trajectory straying into this region non-singular. To see 
this, consider the addition of a shell with a source four- vector given by Eq.(4.1) at i = £o. 
Both a;+ and a;_ will suffer a discontinuity at the shell. The discontinuity (Aa;±) is given 
by integrating Eqs.(3.1a) or (3.1b) across the surface: 

Aio± = -S7rR{eo){aTj). (3.15) 

By a suitable choice of a and j it is always possible to raise or reduce one or the other of 
u± while leaving the other unchanged. In particular, as the arrow on the point Q indicates 
the value of uj- can be reduced in such a way that the trajectory is delivered back to safety 
albeit by flirting dangerously close to the singular point Q. 

What is the physical significance of the points V and 7^'? These are both fixed points 
of UJ+, UJ- and R outside matter: 
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a;' = uj'_^0 R' = 0. 



(3.16) 



As a result, R" = and all higher derivatives vanish at these points. If we exit matter at 
any one of these two points, the spatial geometry degenerates into an semi-infinite cylinder, 
X outside. They are clearly unstable fixed points. Under any small perturbation, 
the vacuum trajectories terminating on either of these points will find themselves either 
returned to the origin, P or consigned to singular oblivion. 

The radius of this cylinder is fixed by the value of the ADM mass. We note that 
Eq.(3.8) implies 



We have sketched the exterior behavior explicitly for a = 2 on fig. (3.1). How sensi- 
tively dependent is this picture on the gauge parameter, a? 

If a is reduced below two, the points V and V' slide out along the R' = diagonal in 
opposite directions reaching infinity at o; = 0.5 — the superspace lightcone value. If a > 2, 
V and V' converge on the R' = diagonal, coinciding asymptotically on the Kji = axis 
as a — > oo, the polar gauge value (discussed in paper I). We see explicitly how polar gauge 
imitates a moment of time symmetry in a very singular way. 

When a = 1, Eq.(3.9) reads R'^ > {RK ji^l)'^ , where the =F correspond respectively to 
\Kfi\ = ±Kfi. These are two cones with apices at the points V and V' given by RKji = ±1, 
R' — 0. Qq is now simply a square, whereas Q is some more complicated figure. 

We note, however, that while the partition of the {uj+, uj-) plane depends qualitatively 
on a, topologically it is identical to the partition illustrated in fig. (3.1) for a = 2. 

In the vacuum region the trajectory is determined by two constants, moo ^ind Cq,, 
which are related to the optical scalars via 




(3.17) 



8m, 



'OO 



(3.18) 



UJ+UJ. 



= 4- 



R 



and 



a;+ — (jj- = ARKji = 



4C, 



a 



(3.19) 
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There are two possible solutions 



2C i?2a-i(i?-2moo) 

^±=:^(±i + \ 1 + ) 



4C„ R^»-^{R-2moo) 



(3.20a) 



and 



4C, / i?2«-i(i?-2moo) 



(3.206) 



a 

From these equations it is easy to analyse the behaviour of a trajectory as it approaches 
a singularity. Let us assume that Cq, > 0. Since we approach a singularity we can assume 
that R is small and positive and that 4i?' = a;+ + a;_ is negative. This means that the 
trajectory is given by Eq. (3.21b). Further, since a > 0.5 the second term in the square 
root is much less than 1 and we can use the Taylor expansion to get 

2mooR'^~^ 
a;_|_ ~- 



-4C, 

UJ- ~ 



(3.21) 



Thus we see again the same structure that we described in Section 2 whereby the 
singularities we get when Ca ^ are stronger than in the MSG configuration. 

4 THE SPHERICALLY SYMMETRIC SHELL WHEN J ^ 

Clearly we cannot solve Eq.(2.4) exactly. Furthermore, if J 7^ 0, even the uniform 
current /density model becomes non-trivial. The only model we will solve exactly is the 
shell. The dynamics of moving shells is a subject which has received extensive study. Our 
focus of interest will, however, be restricted to an examination of the constraints and the 
identification of constraints on the sources avoiding singularities. From one point of view, 
we have already essentially solved the problem in our examination of the exterior solution. 
This is because all of the interesting physics occurs in this exterior region. 

In paper II, we examined the corresponding MSG. A very rich configuration space is 
revealed when we relax J = 0. Let 
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where 



j=jSie-eo), 



(4.1) 



In this form, ctq is the rest mass of the shell. If ctq is real we satisfy the DEC. 

Inside the shell, the space is flat so that R = £. As was the case at a moment of time 
symmetry the material energy M coincides with its Newtonian value 



M = Airall , (4.2) 



and is unaffected by the motion of the shell. 

In any a-gauge, the momentum constraint implies 



,0 £ < £o , , 



In particular. 



Kfl(4+) = 47rj (4.4) 

is independent of a. Kn is finite everywhere so long as R remains bounded from below out- 
side the shell. We can now integrate Eq.(2.1) across £ = Iq io determine the discontinuity 
in R' at the shell 

Ai?' = -AttctIq . (4.5) 

This discontinuity is independent of J. Eq.(4.5) implies that on the outer surface of the 
shell, 

i?'(4+) =1 - 477(74 

_ _M (4.6) 
If (To > 0, i?'(£o+) is bounded above by one. The ADM mass is now given by 

moo = m(4+) = I (l - ^"(4+)) + \^lKl{£o+) • (4.7) 
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We now substitute Eqs.(4.4) and (4.6) into Eq.(4.7). Exploiting Eq.(4.2), we can express 
moo in terms of M, £q, a and j as follows, 

moo = M-—il-v^), (4.8) 
where we introduce the notation 




We note that 

is manifestly positive if the DEC holds. We note that a does not appear in Eq.(4.8). 
Eq.(4.8) generalizes Eq.(7.6) of paper II. The binding energy defined in I, 

/U-2 

-EB = M-m=—{l-v^), (4.80 

is diminished below its Newtonian value by the motion of the shell. It is, however, still 
negative whenever matter satisfies the DEC, consistent with our hopes and allaying our 
fears. If ctq = corresponding to a null shell which saturates the DEC (moving either 
inward or outward), Eb = 0. 

The divergence of outward bound future (past) directed null geodesies at the surface 
of the shell is given by 

(e±)o = I (l - 47r£of7(l T v))) . (4.9) 

We see that (©+)o < and hence that a future horizon must form at some point outside 
the shell whenever 



47Ta£o{l -v)>l. (4.10) 

If V = 1, corresponding to a null outward moving shell, no future horizon can form.* As 
one would expect it is easier to form an horizon when v is negative. 
If AR' < —2 the geometry will be singular. This reads 



This does not, however, mean that the geometry cannot be singular. 
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which is independent of v. 

Let us interpret the exterior solution we examined in sect. 3 in terms of the parameters 
of the sheU model. We note that moo is given by Eq.(4.8) and 

4.11 

If moo is negative, the positive mass theorem tells us that the geometry must be singular. 
We note that moo < in Eq.(4.8) implies 

The potential V{R) is monotonic and unbounded from below. We note that Eq.(4.12) 
implies that R'{io+) < -1 (see Eq.(4.6)).t 

If moo is positive, V{R) possesses a maximum at the point Rc given by Eq.(3.6). 
There are two factors which determine the nature of the exterior geometry. In sect. 3 we 
saw that the condition V{Rc) < 1 can be cast in the form (3.8). Substituting Eq.(4.8) for 
moo and (4.11) for Cq,, we get 

(l - |(1 -«^))" < (2a- 1) (^)"„^y^(-«>. (4.13) 

where y := M/Iq. This relationship clearly depends on the choice of a. When a = 1, the 
condition V{Rc) < 1 is simple: 

M 1 

> 



24 1 + l^^l 

If, in addition, i?'(£o+) < the geometry will be singular. Eq.(4.6) then implies 

^>1. (4.14) 
When the DEC is satisfied all such geometries are singular. 



t It will be a strong J-singularity unless v = which is the only way that C can 
vanish. The simplest model with a non-trivially vanishing C consists of two shells moving 
in opposite directions. 
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We also note that £o < Rc is equivalent to 



M 1 

> 



24 1 + (2a - 1)^2 • 



When a = 1, 



M 

> 



But, when the DEC holds, 



2£o 1 + 



1 1 

> 



1 + ~ l + \v\' 

There are therefore no geometries which simulataneously satisfy V{Rc) > 1 and £0 < Rc 
when the DEC is satisfied.* Thus all geometries with V{Rc) > 1 are non-singular when the 
DEC holds. We illustrate the situation in fig. (4.1). We are now in a position to conclude 
that 

M<^^. (4.15) 
- 1 + |i;| ^ ' 

in any non-singular shell geometry. 

Later, we will discover that the (weaker) bound. 



M < 24 , (4.16) 

holds in a non-singular geometry regardless of J and as a result also of the DEC. We recall 
that the shell saturated this condition when v = (the MSC result). 

It is useful to recall the moment of time symmetry analysis. It was conjectured by 
ADM and subsequently proven by BMOM that [8] (also see ref.[2]) 

M -m> M^/2eo. (4.17) 

The conjecture was motivated by the fact that in Newtonian gravity, the configuration 
that minimized the binding energy for a given total M is the shell. If (4.17) holds, then 
Eq.(4.16) follows by the positivity of the quasilocal mass in any non-singular geometry. It 
is important that the bound (4.16) is saturated by a shell. 

* Equality £0 = Rc gives R' = outside so that the exterior of the shell is a cylinder of 
radius £0. This is singular. 
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When V 7^ 0, Eq.(4.15) is stronger than (4.16). When J ^ 0, we would expect the 
analogue of Eq.(4.17) to imply a bound on M at least as good as (4.16). 

What is this analogue? If we were to take Eq.(4.8') at face value, we would conjecture 

that 

M-m> , (4.18) 

where P = Att J jE^dl = AttJIq. The positivity of m, however, would now imply that 

M(l - P'^/M^) < 2io , (4.19) 

which is considerably weaker than Eq.(4.16). The problem is that in the shell, when v 
is large the only solutions with m simultaneously small and positive are all singular (see 
fig. (4.1). To be more precise, if we substitute the inequality Eq.(4.15) into Eq.(4.8') we 
get 

m > M\v\ . 

Hence m is bounded from below if v is non-vanishing. 

This suggests that we can do better than (4.16). In fact, Eq.(4.15) suggests that 

M + \p\<2e. 
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THE GLOBAL CHARACTERIZATION 
OF TRAPPED SURFACES AND SINGULARITIES 

At a moment of time symmetry, there is a remarkable similarity between the signal 
for the presence of an apparent horizon, R' = and that for the presence of a singularity, 
R = 0. In paper II, this meant that the techniques which were good for analyzing apparent 
horizons were almost always also good for singularities, and the effort required almost 
identical. In general, however, the signal for an apparent horizon will involve the extrinsic 
curvature of the spatial hypersurface and we need to distinguish between future and past 
horizons. It is this feature which complicates the analysis of apparent horizons. It is 
remarkable that the non-triviality of the momentum constraint and its coupling to the 
Hamiltonian constraint does not present a serious obstacle. 

5 SUFFICIENCY 

In paper II, we demonstrated that sufficiency conditions for the presence of trapped 
surfaces and singularities at a moment of time symmetry could be cast in terms of inequal- 
ities of the form, if M > some constant times io, the geometry must contain a trapped 
surface for one constant and a singularity for some other constant. In this section, we 
generalize the inequalities of this form to general initial data. 

As we will see, the natural generalizations of M which arise are the quantities M ±P, 
where 

P:=Att / MR^J , (5.1) 

JO 

is a measure of the total current. We note that when the DEC is satisfied, then 

M±P = Att MR^{p±J) (5.2) 
Jo 

is positive. A few years ago, Bizon, Malec and O Murchadha demonstrated that if 

M-P>^£, (5.3) 

on a maximal slice, assuming only that p > 0, the spatial geometry must contain a future 
trapped surface [5,6]. This generalizes the inequality which is valid at a moment of time 
symmetry. They also showed that the numerical coefficient appearing on the RHS is 
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sharp. They did this by constructing a solution with M — P > (7/6 — e)£ but without any 
trapped surface. This solution notably did not satisfy the DEC. More recently, Malec and 
O Murchadha were able to prove that if the DEC holds, the improved bound, 

M-P>e, (5.4) 

holds [7] . They did this exploiting in a striking way their reformulation of the constraints 
in terms of the optical scalar variables. The inequality (5.4) is particularly impressive 
because it coincides with the MSC result when P = 0. 

Unfortunately, neither of the inequalities (5.3) and (5.4) involves spacetime scalars 
on the LHS so it is not clear what invariant significance they possess. Does a change of 
foliation change these results? To examine this question, in this section we will examine the 
issue within the framework of the one parameter family a foliations. We will inequalities 
similar to (5.3) and (5.4) which are not tied to maximal slicing but are valid in any valid 
a-gauge. 

5.1 Trapped Surfaces: Weak Energy, a - Gauge 

Let us first assume that only the weak energy condition is satisfied but instead of 
considering only a = 2 as BMOM did in their derivation of Eq.(5.3), we will suppose that 
0.5 < a < oo. When a = 2, Eq.(5.3) is satisfied. We will prove that, in general, 

M - P > f{a)£, (5.5) 

where 

f(a):=l + - ^ ' , (5.6) 

reproducing the bound Eq.(5.3) when a = 2. However, the minimum of /(a) is assumed 
when a — 1 where we reproduce Eq.(5.4). Curiously, the gauge providing the best bound 
when we do not assume dominant energy is not maximal slicing. The likely reason for this 
is that in this gauge, = P/R^. 

The original proof by BMOM exploited conformal coordinates. Our approach eschews 
tying ourselves to any particular spatial coordinate. Not only is the end result independent 
of the spatial coordinate, it is clear that the coordinate invariant approach is not only more 
transparent but also more efficient. 
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When Eq.(2.3) holds we can rewrite Eq.(2.1a) in the form 



47rpi?2 + d,{RR') = i (l + (R'f) + ^-^R'Kl . (5.7) 
We integrate from £ = up to the surface value £: 

M + RR' = r+^—^ I dlR^Kl, (5.8) 
where T is defined (as in paper II) by 



1 

r := - / de 

2 Jo 



l + {R'f . (5.9) 



We now eliminate R' in the surface term in favor of the optical scalar uj+ and Kr using 
the defining relation (2.6a). The vanishing of a;+ signals that the geometry possesses an 
apparent horizon. Let us assume that the surface is not trapped so that a;+ > 0. 

To eliminate the Kr dependence on the boundary which comes along with the replace- 
ment of R' by u>+ in Eq.(2.6a), we note that we can integrate the momentum constraint, 
Eq. (1.3.2) to obtain* 

R^Kr^ P + {l-a) / MRR'Kr. (5.10) 

JO 

Substituting Eq.(2.6a) and (5.10) into (5.8) we now obtain 



M 



1 

- P + 2a;+ = r + - (1 - 2a) / dlR^K\ - {1 - a) / dlRR'KR . (5.11) 
2 Jo Jo 



When a > 0.5, the second term on the RHS is manifestly negative. As such we could 
discard it to cast (5.11) as an inequality. However, it is clear that we can do better by first 
completing the square in the sum of the second and third terms before discarding: 



1 

-{l-2a) / dm^K\ - {I- a) j diRR'KR 

1, - [^{uL . l^U')\ 1(1=^ /' ^U'^ (5.12) 



2 2a 



The privileged role of the gauge with a = 1 is evident. 
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On the last line, we have used the fact that when p is positive and a > 0.5, R'^ < 1. In 
addition, under these conditions, we obtain the same upper bound on F, 

T<e. (5.13) 
as we obtained at a MSG. We conclude that if the spherical surface is not trapped, then 

l(l-a)2" 



M -P< 



1 + 



i. (5.14) 



2 2a-l 

Thus if the surface is trapped or the interior contains a trapped surface Eq.(5.5) holds. 
Let us now examine some extreme cases: 

We note that / diverges as we approach the minisuperspace lightcone, a — > 0.5 and 
a ^ oo. While it is tempting to interpret this as a signal of the breakdown of the gauge 
on the lightcone, it is also clear that the discarded negative term blows up at these two 
values. 

Let us consider the two extreme distributions saturating the DEC everywhere, P = 
±M which are respectively the cases of a radially outward and a radially inward moving 
null fluid. In the former case, Eq.(5.14) becomes a vacuous statement — even though we 
do expect it to be more difiicult (if not impossible) to form an apparent horizon. In the 
later case, we have that if 2M > f{a)£, the geometry will possess a trapped surface. It 
is twice as easy to form an apparent horizon with an inflowing null fluid as it is with a 
stationary fluid. 

It is also possible to tighten the sufficiency condition in the same way we did for MSCs 
when p' < if, in addition, J has a fixed sign. 

We note that in same way we did when Kab = 0, when p' < we can replace Eq.(5.13) 

by 

^ <i-^-^ [ MR^JKr- [ MR^Kl. (5.15) 
o R Jo Jo 

If J is positive (negative) everywhere, the inequality still holds when the third term on the 
RHS is dropped. The (negative) last term on the RHS can now be added to the (negative) 
term of the same form in Eq.(5.11) before the completion of the square. We get 



4M 



P < 



2 



27 



As before, this is minimized when a = 1 and when P = again reproduces the result at a 
moment of time symmetry. We note that both the LHS and the RHS have been improved. 
When a = 0.5, unhke Eq.(5.14) the RHS of Eq.(5.16) does not diverge. From one point of 
view, Eq.(5.18) is not very satisfactory — we have broken the symmetry between J and 
p. However, on the other hand this asymmetry permits us to write down a non-vacuous 
sufficiency condition when the the spatially averaged DEC is saturated with P = M. 
Whereas Eq.(5.14) is vacuous under these conditions, Eq.(5.18) provides the non-trivial 
statement: suppose p' <0 and the motion of matter is outward and null, then if 



M > 3 



1 + 



2(2a + 1) 



(5.17) 



the spatial geometry will possess an apparent horizon. When a = 1, this value of M is 
three times larger than that of a corresponding stationary distribution of matter. 

5.2 Trapped Surfaces: Dominant Energy, a - Gauge 

When the DEC holds, our experience suggests that the appropriate variables are the 
optical scalars. The optical scalar which marks the presence of a future trapped surface 
is u>+. Remarkably, only the constraint (3.1a) determining the spatial derivative of a;+ 
will play a role in the determination of the inequality. Let us first recast Eq.(3.1a) as an 
equation for the spatial derivative of Riu+: 

{Bw+y = -8nR^{p -J) + ^ [2u;+u;- + 4 + 4:RKu+ - loI] . (5.18) 
This equation can be integrated up to give 

Ruj+ = -2(M -P) + 2T+ , (5.19) 

where 

1 

r+ := - / d£ [2u}+u}_ + 4 + 4:RKu}+ - a; J] (5.20) 

8 Jo 

are the natural optical scalar generalizations of F. In particular, when Kab = 0, r_|_ = F. 
In general, in the gauge (2.3), 



RK ={2- a)RKR = (2 - Q;)-(a;+ - a;_) , 
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(5.21) 



so that 



1 r 

r+ = - / M[A + auj+uj- + (1 - cx)u\] . (5.22) 
8 Jo 

X 

Let us examine two special cases. The case examined by MOM was a = 2. Now 

1 

r+ = - / d£ [4 + 2a;+a;_ - u\] . (5.23) 
We note that < (a;+ — a;_)^ imphes 

2a;+a;_ - a; J < a;^ . (5.24) 

We now exploit the inequality (2.19) to obtain r_|_ < I. If the surface is not future trapped 
then a;+ > in Eq.(5.19) reproducing Eq.(5.4). If cu = 1, then* 



1 



+ co+uj-]<e, (5.25) 



using Eq.(2.22), which is the same as the a = 2 value. The a = 1 bound does not improve 
even though the energy condition is more stringent. This suggests that this bound is sharp. 

What if a is not one of these two special values? In general, the bound on both cv+ 
and U- depends on K. This makes it less obvious how to bound r_(_ for any a other than 
these two values. What we can do is bootstrap on Eq.(5.21) to turn this into a bound 
which is independent of K. Let kr = Ma,x\RKji\. We note that Eq.(5.21) implies 

«fl < ^(Max|a;+| +Max|a;_|) . (5.26) 

We know from Eqs.(6.4a & b) in paper I that Max|a;_|_|) = Max|a;_|) < k + -^/[/tpT^. 
Hence 

<l(AC + (|ACp+4)V2) 

1 , , (5.26') 

= - (|2 - a\KR + (|2 - + 4)^/^ J . 

It is straightforward to invert Eq.(5.26') to obtain the bound 

KR < , . (5.27) 



a\ 



Note that in this gauge, r+ = r_ . 
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This bound unfortunately is valid only for 1 < o: < 3. Eq.(5.21), in turn, implies the 
bounds on u>+ and a;_ , 



We note that when a = 2 we reproduce the bounds, Eq.(2.19). What is remarkable is that 
these bounds are independent of J. Note that our knowledge of the bound on R'"^ does 
not help (nor should it be expected to help) to improve these inequalities. 

We now return to Eq.(5.22). We introduce a parameter b, and we rewrite 



+ 
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a + (1 — a)b^LO+u;- + (1 — Q:)(a;^ — bu+LO-) 



(5.29) 



We now complete the square on the last term and use an obvious modification of Eq.(5.24) 
to obtain the inequality 



r+< 



d£ 



4 + (cK + (1 - a)b)iO-i-iO- + («-!)- 



-U)_ 



(5.30) 



We can now exploit Eqs.(2.22) and (5.28) to obtain. 



r+< 



1 + q;+ (l-a)6 + 
l + a+{l-a)b+^ 



i, l<a<2 
2 < a < 3. 



1 

4 



(5.31) 



The idea is to find the b which minimizes the RHS. For fixed a, the RHS is minimized 
when 



2{a - 1) 
2(3 -a) 



1 < a < 2 

2 < a < 3. 



(5.32) 



The corresponding values of r+ are 



< 



itt(3-a)£, l<a<2 
lla"^ -3a + 4]i, 2 < a < 3 . 



(5.33) 



These bounds reproduce the optimal values obtained when a = 1 and a = 2. In the 
neighborhood of a = 2, this bound is an improvement over the bound (5.5) which does not 
assume the DEC. Even though the bound (5.27) diverges both as a — > 1+ and as a — > 3_, 
the bounds (5.14) at these limit points nonetheless are finite, and in the former case as we 
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have just seen even coincides with its optimal value there. It is not clear how to extend 
this technique outside the range a = 2 ± 1. 



5.3 Singularities 

It is not obvious how to import the DEC into the statement of a sufficiency condi- 
tion for singularities. What we have is an obvious generalization of the moment of time 
symmetry result: We recall that, in general, 



We proceed exactly as for a moment of time symmetry. F is always bounded by one 
whenever p is positive. Furthermore, i?' < 1 so that R{t) < I everywhere on a non- 
singular geometry and B! > —1. The surface term is therefore bounded from below by —£. 
Finally, the second term is negative whenever a > 0.5 — The Kab dependence is trivially 
handled. Thus we get 



independent of the value of a > 0.5 which is exactly the result at a moment of time 
symmetry. 

As at a moment of time symmetry, if we place constraints on the sources it is possible 
to tighten the inequality. We note that when p' < and J is positive (or negative) 
everywhere, then Eq.(5.15) can be truncated even more brutally, F < £ — M/3. We get 



Unlike the moment of time symmetry discussion we cannot claim that this represents a 
universal bound when p' < and J is positive (negative). The reason is that the geometry 
can still turn singular if J is large enough. 

We note that there is no obvious way of introducing P into either Eq.(5.35) or 
Eq.(5.36). The singularity condition is not symmetrical in M and P. 

6 NECESSITY 

We noted in paper II, in our examination of a moment of time symmetry, that the 
necessary conditions we were able to formulate with respect to M and £ were extremely 




(5.34) 



M <2£, 



(5.35) 




(5.36) 
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weak. If J ^ even these conditions appear to be beyond our reach. What one can do is 
provide generahzations of the necessary conditions which where formulated with respect 
to the variables, pMax and io. These inequalities assumed the form 



PMax^^ < constant. (6-1) 

Typically, we would expect | JMax| and a to enter into this description. We would expect 
that by appealing to the DEC the inequalities should simplify. Crucial to the derivation 
of Eq.(6.1) are two simple Sobolev inequalities of the form 

rh nil 

S / diR^ < / diR'^, (6.2) 
Jo Jo 

where S depends on the boundary conditions satisfied by R. In general R{0) = 0. At 
the first trapped surface, R'{£i) = and S = ti"^ At a singularity, -R(^i) = and 
S = n"^ / (\. At a singularity, we found that R tends to zero like i? ~ (£ — so that R' 

diverges like (£ — ^i)"^/"^. Even though R' diverges so that the integrand on the RHS of 
Eq.(5.9) diverges, the integral itself remains finite. When J 7^ 0, however, R diverges more 
strongly, R ~ (£-£1)^/^+" (see Eq.(2.14)) so that R' ~ - £1)-"/^+". Thus the integral 
on the RHS of Eq.(5.9) will only exist if a < 1 — outside the range found to provide the 
best sufficiency results in Sect.5.1. Thus, whereas we found that we could optimize the 
inequalities of necessity at a moment of time symmetry by weighting R'^ by an appropriate 
power of i?, a non-trivial weighting will be essential when J 7^ at least in the case of 
singularities. 

A Bound on Kn 

To form a necessary condition for singularities it is important to possess some control 
over Kji in a manner which does not require the geometry to be regular. In particular, we 
cannot exploit Eq.(5.27) which is only true in regular initial data. It is, however, simple 
to obtain a bound on Kn by | JmoxI without making any assumptions about the regularity 
of the geometry. We have that 

In general. 
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l^«l<^^^/°^^'^"- (6-4) 

This is the result we will exploit below. There are some interesting related inequalities. 
Suppose that the DEC holds, and a = 1. We obtain 

M 

\Kn\ < ^ . (6.5) 

This inequality in turn implies that the proper spatial average of \Kji\ is bounded by the 
product, pMax^i: 

< li^fll > <47r pMax^l, (6.7) 

a pretty result, even if we have not found an application for it. 
A Bound on R' 

We will also require a bound on R' which does not require the geometry to be globally 
regular. To obtain this bound, we note that in any a-gauge, the Hamiltonian constraint 
Eq.(2.1a) reads 

RR" = ^(1 - i?'^) + ^i?^(l - 2a)Kl - AttR^p . 

At the origin, we have R' — 1. At a singularity we have R' < whereas at a globally 
regular solution we have i?' — 1 at infinity. If R' has an interior maximum then R" 
vanishes there. Hence at that point we have 

1-R'^ = R^{2a - 1)K\ + SttR^p . 

Thus, if we have a standard a-slice, i.e., a > 0.5 and if the source satisfies the weak energy 
condition, p > we must have 

1 - < , 

and therefore at the maximum of R' we must have R' < 1. Therefore this is a global bound 
independent of whether the slice is regular or not. 

6.1 Singularities 
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The most naive generalization of Eq.(6.1) would be an inequality treating pMax and 
I -/Max I symmetrically, of the form: if 

(pMax + I AlaxD^o < C , (6.8) 

for some constant c, the geometry is regular. However, our experience examining the 
approach to singularities suggests that this is too optimistic. The natural inequality we 
obtain involves not J but its square, assuming the form: if 

(pMax^O + Cl(| JMa^Ko)^) < ^2 , (6.9) 

where Ci and C2 are two constants, the geometry is regular. Even if matter satisfies the 
DEC, once we foliate extrinsically the symmetry is broken. The value of J plays a more 
significant role than the value of p. This is consistent with our findings in Sect. 2 in our 
examination of the generic behavior of the metric in the neighborhood of a singularity 
in an cu-foliation of spacetime. The optical scalar variables suggest that a more judicious 
gauge involving some mix of intrinsic and extrinsic variables might restore the symmetry 
between p and J we have broken with the a-parametrized gauges. 
We note that Eq.(2.1) implies 

^(1 + R'^) = {RR'Y + AirpR^ + ^(2a - 1)R^K\ . (6.10) 

The last term is manifestly positive. Suppose that the geometry is singular oX I = l\. We 
cannot simply integrate Eq.(6.10) and discard the boundary term. First of all, R!^ is not 
integrable on the interval [0, £i] and, secondly, the surface term does not vanish at 
the singularity unless a < 1. 

What we need to do is multiply Eq.(6.10) by some (positive) power of R before inte- 
gration. The relevant power of R will generally depend on the value of oi. To restore the 
divergence appearing in Eq.(6.10) we need to perform an integration by parts. We now 
integrate up to i\. 



1 rt^ 
'2 Jo 



(MR^il + (2a + l)i?'2) = + 47r / MpR^^" + -(2a - 1) / R^+^'K 



(6.11) 

To discard the boundary term, we require R^^°'R' to vanish at the singularity. This implies 
that 
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a>a-l. (6.12) 

This choice of a simultaneously bounds the integral over R"^R'^. 

We also will need to place a bound on the last term on the RHS of Eq.(6.11). We 
exploit Eq.(6.4) to bound Kr. The problem is that this bound involves the positive power 
of i?, in the denominator which is difficult to control. We obtain the bound, 

^ ' diR^+^'Kl < {AirfJM^j^ ' cZ£i?"-2" , (6.13) 

on the term quadratic in Kji. If the weighting term is chosen such that 

a>2a, (6.14) 

the denominator problem is solved. Fortunately, such values are consistent with Eq.(6.12) 
for all physically acceptable values of a. The RHS of Eq.(6.11) is clearly simplest when 

a^2a. (6.14') 

This is the value we will henceforth adopt for a. The expression is still not very useful as 
it stands. A remarkable fact, however, is that we can bound it by an integral over R'^^^'^°'\ 
To understand why this is important, note that the integral over i?'^ appearing on the 
LHS of Eq.(6.11) can be cast in the form 

/ MR^'^R'^^-, — r / ci£ (i2"+^)'2 . (6.15) 

Jo (« + 1) Jo 

The Sobolev inequality can be exploited to place a bound on the integral over the function 

^2(l+a). 



Jo 



1 /-^i 

^^2(a+i)<^ — diiR^'+^y^, (6.16) 
{a+ ly Jo 



where the constant Sq = 7r^/^f is the Sobolev constant which is relevant for functions 
which vanish at both £ = and £ = £i. 

We now prove the existence of a bound of the form 

rii /re \^ i>ei 

J Mi J MR^+'^j <CJ c/£i?2(«+i), (6.17) 
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for some appropriate constant C. A crude bound is provided by the positivity of the 
covariance for any power i?"': (Holder Inequahty), 



2n 



(6.18) 



which imphes 



(6.19) 



so that 



rii / ri 

J diij MR^ 



(6.20) 



We can, however, do better. Let 



Now G'(O) =0 and G'{li) 
the appropriate constant 



G{t) := / d^R'' . (6.21) 

JO 

= 0, for all n > 0. We apply the Sobolev inequality to G with 



so that 



(6.22) 



(6.23) 



This is better by a factor of 7r^/8 than the estimate (6.20). The end result is the bound 

(6.24) 



JQ Jo 



on the third term on the RHS of Eq.(6.11). We can now write 



1< 2 



47rpMax + 32(2a - 1) J^^.^? - ( - 



TT \ 1 + 4a 



£J 2(l + a)2 



^^2(l+a) / / ^^2a 



(6.25) 
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In paper II, we proved that the ratio of integrals appearing on the RHS can be bounded 
as follows (Eq.(II 6.3.16)) (a = 2a) 



J^R'^de ~ 3 + a ' ' 



(6.26) 



which implies 



1 3 + 2a 1 + 4a 



+ 



2 1 + 2q; 2(1 + q;)2 



tt" < ^npMaJl + 32(2« - l)Ji^^jt 



(6.27) 



If a = 1, 
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1 + 



37r2 



'2\2 



< H (<^Max-^i) 



(6.28) 



In Eq.(6.27), it does not make much sense to claim that one value of a provides a 
better bound than another value. Not only is the inequality invalidated if the quadratic 
in JMax is dropped, Jnax plays a more decisive role than pMax in the inequality (6.27), 
appearing as it does through its square in contrast to p which appears linearly. The MS 
inequality does not generalize in the obvious linear way. 

If the DEC, the inequality simplifies. For a = 1 we obtain 



- H < 32 PMax-^? H 

6 8 ~ y ^ 16/ 8 



(6.29) 



or 



5 3 „ TT 
1 TT^ 

3 2 2 



< PMax-^l 



(6.30) 



The LHS ~ 5/16, which is approximately half as good as the moment of time symmetry 
result. 



6.2 Apparent Horizons 

We note that Eq.(2.1) implies 

-(1 + R''^) = {Ruj+ - R^KrY + 47rpi?2 + l(2a - 1)R^K'^ . (6.31) 

Again the third term on the RHS is manifestly positive. We can integrate Eq.(6.31) up to 
the first future horizon 
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r = -RKR\e^ + 47r / ' MR^p + ^(2a - 1) /" ' MR^K\ , (6.32) 
Jo 2 Jo 

where F is given by Eq.(5.9). We note that Eq.(6.3) places a bound on Kr in the surface 

term. Thus 

r < 47r(pMax + \JMa.\) J^^ + ^(2« " 1) / ^ . (6.33) 

A linear term in Jm-ax appears in the apparent horizon inequality condition which is not 
present in the singularity inequality. This is a reflection of the different boundary conditions 
enforced there. 

We can exploit a Sobolev inequality to place a bound on the integral over the interval 
(0,£i) of the quadratic R^ by the same integral over the quadratic, R'^: 

S MR^ < MR'^. (6.34) 
Jo Jo 

The relevant boundary conditions are 

R' + RKr = (6.35) 
at £ = £i . The inequality is saturated by the trigonometric function, 

R{t) = sin(7£) , 

which also determines the optimal value of 5" = 7^. The boundary condition, (6.35) 
determines 7 to be the lowest solution of the transcendental equation. 



tan7£i = --^. (6.36) 



We note that 



7 < ^ (6.37) 

if Kr is negative with 7 — > n/2£i as Kr which is the moment of time symmetry 
bound and 7 Tr/ii as Kr — > +00. 

When we attempt to bound the right hand side we run into the same problem we 
faced before with the last term. In addition, however, we must contend with the surface 
term. 
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The same weighting we found worked before works again. To restore the divergence 
appearing in Eq.(6.31) we need to perform an integration by parts. We integrate up to £i: 



1 /-^i 



2 



'0 



+ 47r d£pR''^^+°'^ + J(2« - 1) r . 
Jo 2 Jq 

(6.38) 

We now exploit Eq.(6.3) to bound the Kr and Kr term as foUows For the former, 

^2(i+a)^^ < AirR^+'^Jua^ / deR^+'' . (6.39) 
^1 Jo 

The weighting process has broken the symmetry under interchange of p and J of the hnear 
terms on the RHS of Eq.(6.38). For the term quadratic in Kr, we again have ((6.12) with 
a = 2a) 



t R^Kl < (47r)2j^^ ( fd£R^ 
Jo Jo \Jo ^ 

We again require a bound on the last term by an integral over This time, however, 

the Sobolev constant is that which is relevant for functions which vanish at £ = but satisfy 
Eq.(6.35) at £ = £i, i.e., = 7^, where 7 is given by (6.36) and (6.37). 

The crude bound we derived before, (6.20), is expected to work better this time. As 
before, however, we can do better. This time we let 

H{£) := dm"" I ' MR"" . (6.40) 

Now -ff(O) =0 and H{£i) = 1 for all n. We apply the Sobolev inequality to H with the 
appropriate constant 

J d£H{£f<(=^j J d£R^^/(^J d£R''j . (6.41) 

so that Eq.(6.23) holds exactly as before and we again obtain the bound (6.24) for the 
integral over Kr. We can now write 
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1< 2 



47rpMa. + 32(2a-l)j2,^,£2_^2 l + 4a; 



2(l + a)2 

+87rJMax^i^" 



(6.42) 



/o ' Jo 

An upper bound on 7 in Eq.(6.42) is provided by its — > 00 limit, i.e., 7r/£i and the 
lower limit is zero. We can again exploit (6.26) to bound the ratio of the integrals in the 
first term of (6.42). In the second term it is clear that 



R 



l+a 



MR 



2a 



is bounded by £f if a < 1. Therefore a necessary condition for the appearance of a trapped 
surface at some proper radius li is that 



47rpMax^? + 4^^^^Max^? + 32(2a - l)Jl,,j\ > i^±|^ 



(6.43) 



CONCLUSIONS 

This paper concludes a series of three papers on the identification of the configuration 
space in spherically symmetric general relativity. We have attempted to provide a coherent 
synthesis of two very different ways of looking at the constraints, one in terms of the tradi- 
tional metric variables, the other in terms of the optical scalar variables. Which description 
is appropriate depends very much on the details of the problem under consideration. 

A very satisfying representation has emerged of regular closed solutions as closed 
bounded trajectories on the {ujj^,uJ-) plane. In this representation, R plays a secondary 
role. We have performed the analysis explicitly in vacuum. We will show elsewhere that 
this plane also provides a very profitable representation of ct-slicings of the Schwarzschild 
spacetime [9]. 

We have presented a variety of necessary and sufficient conditions for the presence 
of apparent horizons and singularities in the initial data. This paper is necessarily more 
open-ended than either paper I or paper II. It is clear that some of the Sobolev inequalities 
exploited in sect. 6 can be sharpened. Indeed, the professional will consider our approach to 
functional analysis extremely heuristic. As physicists, however, we are more interested in 
the fact that such bounds can be established than in squeezing them for better constants. 
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Where does one go from here? The obvious challenge is to generalize this work to 
non-spherically symmetric geometries. One needs to bear in mind, however, that our 
ability to describe the configuration space in considerable detail has relied on features of 
the spherically symmetric problem which we know do not admit generalizations. 

There is still, however, much that needs to be done before we can claim to understand 
spherical symmetry. 

We need first of all to examine the classical evolution. Write down the Einstein equa- 
tions with respect to the optical scalar variables. Can we cast the theory in Hamiltonian 
form? If the value of these variables in the analysis of the constraints is anything to go by, 
one has every reason to expect that they will throw light on the solution of the dynam- 
ical Einstein equations, both analytically and numerically. Indeed Rendall has recently 
exploited these variables to extablish a global existence result [10]. 

A physically interesting question that is extremely relevant is the identification of 
initial data that potentially might develop apparent horizons. In principle it should be 
possible to do this exploiting in addition to the constraints the dynamical Einstein equa- 
tions evaluated on the initial hypersurface. These equations involve the pressure of matter 
though some equation of state. The scenario which is most susceptible to collapse is 
pressureless matter. We should be able to exploit this condition to formulate neceasary 
conditions along the lines developed in sect. 6. At the other extreme, a stiff equation of 
state would inhibit collapse. Thus such a scenario might provide a sufficient condition. A 
successful analysis of this nature has the promise of putting an analytical handle on the 
physics hinted at in Choptuik's numerical simulations of the collapse of a massless scalar 
field [11]. 

Finally, the bounds on the optical scalars are certain to have profound implications 
for the canonical quantization of this model for gravity [12]. We hope to examine this 
problem in a subsequent publication. 
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Figure Captions 



fig.(3.1) Non-singular exterior vacuum solutions on the {iv-^-,uj-) plane for a = 2. All 
non-singular exterior trajectories lie within the 'wedge' shaped region, Qq, bounded by the 
arc segments, CQ and QD. 

fig. (4.1) M/2£q vs. V in the shell model. All non-singular geometries lie below the curve 
V{Rc) = 1. 
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